Tunable Fano resonance in a parallelly coupled diatomic molecular transistor 
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We investigate electron transport through a diatomic molecule parallelly coupled to infinite source 
and drain contacts. We utilize a model Hamiltonian involving a Hubbard term in which the contacts 
are modeled using recently developed complex source and sink potentials. The zero bias transmission 
spectrum for a symmetrically coupled system as a function of the Fermi energy acquires a Fano 
lineshape as the Hubbard interaction is turned on. For large values of U, the Fano lineshape broadens 
and shifts to higher energy values disappearing eventually. Meanwhile, the Breit-Wigner resonance 
located at the bonding resonance in the noninteracting limit survives but its position is shifted twice 
the coupling between the atoms in the molecule in the infinite U limit and its linewidth is reduced 
to half. We attribute this behaviour to the unavailability of one of the transmission channels due to 
Coulomb blockade. 
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Molecular electronic ;! 1 ' 2 ' 3 ' 4 ' 5 -' was put forward as a novel paradigm for the construction of future electronic devices 
almost three decades agoP However, its feasibility has been much debated since it seemed to offer no clear alternative 
to the conventional bits for the transmission of information. The advent of spintronics^', where spin and charge 
transport are investigated in mesoscopic systems, provided a satisfying resolution to this controversy by proposing to 
utilize the spin degrees of freedom of the electron as bits. These developments culminated in the emergen ce of the 
field of molecular spintronics which aims to incorporate spin dependent transport into molecular electronics! 9 ' 10 ! 

Effective one-electron theories based on a combination of non-equilibrium Green's function tech nique and density 
functional theory have dominated the theoretical analysis of molec ular e lectronic devices (MED's) 11 ' 12 ' 13 ' 14 ' 15 ' 16 ' 17 ! 
However, these models fail to account for electron correlation effect s^ 7 -^ since the choice of an exchange-correlation 
functional which can capture the self-energy of the electron adequately is unknown. This is reflected in large dis- 
crepancies between theoretical predictions and experimental results. Therefore, it is of crucial interest to be able to 
develop theoretical methods that can ta ckle t his problem. 

Hedin's celebrated GW approximation^' 2 ^' has been invo ked rec ently in conjunction with a model Hamiltonian to 
overcome the self-energy defi cienc ies of one-electron theories! 21 ' 22 ' 23 ! For the same purpose, we implemented the source- 
sink potential method (SSP)!M2f>l m a moc j e l Hamiltonian involving a Hubbard term to probe the strongly correlated 
regime in molecular transistors systematically. In the case of a diatomic molecule serially coupled to ferromagnetic 
contacts, we found that the electron transport is suppressed in the U — ► oo limit while the energy gap between the 
bonding and antibonding resonances is reduced due to a coupling between the molecular singlet and triplet states for 
small C/P 

In this paper, we extend our analysis and analyze a diatomic molecular transistor parallelly coupled to infinite 
ferromagnetic contacts. Starting from the known results in noninteracting limits, we explore the transmission of 
the system in infinitesimal bias as a function of the Fermi energy of the contacts for various values of the intrasite 
interaction strength U. Our results indicate that turning on the interaction results in appearance of a Fano resonance 
when the molecule is symmetrically coupled to contacts. Increasing the interaction strength produces several effects. 
The Breit-Wigner resonance located at the bonding level in the noninteracting limit starts moving towards the 
antibonding level. During this shift, it transfers some of its linewidth to the Fano resonance which also starts moving 
towards higher energies. Hence, sweeping the interaction strength provides tunability for the Fano resonance. Even 
though tunable Fano resonances have been reported for other systems ranging from metallic nanostructures consisting 
of a disk inside a thin ring^ to localized Bose-Einstein condensates in one-dimensional optical lattices^, this paper 
constitutes the first prediction of a tunable Fano resonance in a strongly correlated molecular transistor. In the U — > oo 
limit, the Fano resonance disappears to high energies whereas the Breit-Wigner resonance survives by settling into 
the antibonding position of the noninteracting system with its linewidth reduced to half. 

The source-sink potential method relies on a Bloch wave impinging on the contact molecule interface from the 
left contact. It is partially reflected and partially transmitted at the junction. The reflection coefficient is found 
by solving the Schrodinger equation for the relevant Hamiltonian and the zero bias transmission is obtained from 
the Landauer-Biittiker formula. In this approach, the semi-infinite contacts are described with a Hiickel model and 
we assume that they are ferromagnetic with parallel configurations, i.e. they only contain spin-up electrons. This 
requirement can be realized experimentallyPESl 

The atoms in the semi-infinite contacts are coupled to their nearest neighbour with strength f3c- The first and last 
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atoms in the right and left contacts are each coupled to the molecular atoms through Pcm ■ The Hamiltonian is given 

by 



where 



H = hi + h i + V\ 



n 



(i) 



£ L n LT + T, R n m + Pcmc^c^ + H.c. + 
PcM^c b] + H.c. + l3 CM c^ m c a] + H.c. + 
PcMC^CbX + H.c. + p M c\^ + H.c. 



(2) 



and 



PcM^ Li c ai + H.c. + (3cMC i Ll c bl + H.c. - 
PcMc\ x c ai + H.c. + Pcmc^^ + H.c. 
PMc\ l c bi + H.c. 



n 



(3) 
(4) 



The operators c^ aa (c aa ) with a=a, b and ct CT (cp a ) with /3=L,R create(destroy) an electron in the molecular sites a 
and b and contacts L and R respectively. The source £l = Pc L 



in the Hamiltonian, where q = arccos ^^§^ J, 



and sink "En = Pee 19 potentials ^ appearing 
describe the infinite parts of single-electron contacts rigorously.^ 
Our finite dimensional Hilbcrt space contains the isolated molecule configurations where the spin up and down 
electrons are confined within the atoms of the molecule. This results in four states. Additionally, we include the 
configurations in which the up-spin electron from atom a or atom b is excited into the contacts. This gives us another 
four states. A detailed description of how we obtain these states has been discussed elsewhere. 2(5 This model is 
presumably the simplest realization of a molecule interacting with ferromagnetic contacts. 
The Hamiltonian in Eq. [I] can then be represented in the space of the configurations as 
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The procedure to obtain the transmission probability is straightforward once we are equipped with the appropriate 
Hamiltonian. It has been described in detail previously. Here we just summarize the method briefly. The transmission 
probability T(e) describes the probability of an incoming electron with Fermi energy e to penetrate the molecule. 
It is related to the reflection coefficient r(e), which is buried within the potential El inside the Hamiltonian, by 
T(e) = l-|r(e)| a . 

The Schrodinger equation involving the Hamiltonian given above becomes 



H(e,r)^ = (E 



N-l 



(5) 



where E N ~ 1 is the energy of the N — 1-electron system resulting from removing the up-spin electron! 29 * 30 ! In this 
paper, we are concerned about the evolution of the molecular ground state, thus we set E 1 ^" 1 — Pm where Pm is the 
energy of the bonding orbital in the isolated molecule. 
The corresponding eigenvalue equation turns out to be 



Det(H(e,r)-(p M + e)I) = 0. 



(6) 
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FIG. 1: The transmission probability T(e) for a non-interacting system as a function of Fermi energy e. Red(solid) curve 
corresponds to (3c = — 10, /3a/ = (3cm = —2 whereas blue(dashed) curve is obtained for (3c = —10, (3m = —2 and /3cAf = — 1. 



Eq. [6] determines the reflection coefficient r in terms of the variable e for a given /3m- We solve this secular equation 
for successive values of the Fermi energy e to obtain the entire transmission spectrum. 

We can now embark on our analysis of the transmission spectrum. We start with the noninteracting limit. This 
will help us to verify the validity of our approach and elucidate the implications of the choice of the parameters since 
the transmission results in this limit have already been reported.^ Fig. [TJ depicts the transmission spectrum for a 
noninteracting system, i.e. U=0, for two different values of the coupling to the contacts (3cm- We consider a system 
in which the atoms in the contacts are tightly coupled compared to the intermolecular and molecule-atom coupling 
therefore we pick (3c much larger than the other two couplings. This choice has been motivated by the fact that it 
enables us to analyze a wider transmission spectrum. 
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FIG. 2: The transmission probability T(e) for an interacting system with ?7=0.5 as a function of Fermi energy e. Red(solid) 
curve corresponds to (3c = —10, /3m = /3cm = —2 whereas blue(dashed) curve is obtained for (3c = —10, /3m ~ —1 and 
(3cm ~ —2. 

A glance at this figure reveals that the spectrum consists of nothing but a Breit-Wigner resonance located at the 
bonding energy of the isolated molecule with linewidth equal to 2 | (3cm I- Therefore, our results are in full agreement 
with the previously reported data.^The disappearance of the Fano resonance at the antibonding energy is a result of 
the emergence of slow transitions between the antibonding state and the contacts. Moreover, the investigation of the 
local density of states in each atom within the molecule shed light on this peculiarity by showing that the lorentzian 
Breit-Wigner resonance located at the antibonding resonance becomes a dirac delta function at symmetric coupling. 
Consequently, the antibonding state gets completely localized inside the moleculeP^ 

After confirming the earlier results in the noninteracting limit with our model, we are now ready to start exploring 
the strongly correlated regime. For this purpose, we gradually turn on the intraatomic Coulomb repulsion within 
the molecule. Fig. [2] shows the transmission spectrum for two configurations with different interatomic couplings in 
the weakly correlated regime. The noteworthy feature here is the appearance of a sharp Fano resonance alongside 
with the Breit-Wigner resonance. Fano resonance is an exquisite interference effect arising from the existence of two 
transport pathways, a resonant one and a nonresonant one.^l As a result, it is quite sensitive to the phase coherence 
of the system. The Fano line shape is given by 



(7) 
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where T is the linewidth of the resonance, Gd is the nonresonant conductance and q is the Fano parameter. It is 
defined as the ratio between the resonant tunneling probability and the direct nonresonant tunneling probability. The 
q — > oo limit yields the Breit-Wigner resonance. 

Our transmission results bear close resemblance to the transport spectroscopy of double quantum dot Aharonov- 
Bohm (AB) interferometers. An Aharonov-Bohm interferometer is in essence a double-slit experiment in which an 
electron tunneling from the left contact to the right contact is split into two waves A magnetic flux $ piercing 
the area enclosed by the two paths induces a change in the relative phase of the amplitudes of the two waves. The 
resulting transmission spectrum has been shown to exhibit a Breit-Wigner and a Fano resonance in the noninteracting 
limit.^ The most striking part of our result is that we can achieve the same transmission spectrum in symmetric 
coupling even in zero magnetic field. The Hubbard interaction inside the molecule is responsible for the decoherence 
in our case instead of the applied external magnetic field in AB interferometers. More recently, taking into account the 
spin-orbit coupling in a noninteracting Rashba quantum dot coupled to ferromagnetic contacts has also been shown 
to result in a transmission spectrum composed of a convolution of a Fano lineshape and A Breit-Wigner lineshape.^ 

One particular intricacy that deserves further attention in Fig. [2] is the separation between the Fano and Breit- 
Wigner resonances. It is clear from this figure that the location of the Fano resonance is pushed to higher energies as 
the coupling between the atoms inside the molecule increases. In fact, when | [3m | exceeds twice | (3cm |, the Fano 
resonance is no longer visible in the transmission spectrum even for very small U values. This makes perfect sense 
intuitively because these conditions mean that the tunneling rate of the down-spin electron between the molecular 
atoms is much higher than the tunneling rate of the up-spin electrons from the contacts to the molecule or vice 
versa. Therefore, the existence of two distinct pathways for the up-spin electron gradually disappears leading to the 
elimination of the Fano resonance from the transmission spectrum. 




FIG. 3: Left panel shows the transmission probability T(e) for an interacting system with U —4 as a function of the Fermi 
energy e. Red(solid) curve corresponds to (3c = —10, Pm = (3cm = —2 whereas blue(dashed) curve is obtained for f3c = —10, 
Pm — —1 and Pcm = — 2. The transmission spectrum in the right panel is obtained by using the same couplings in left panel 
for U=S. 

We keep increasing the strength of the electron-electron interaction inside the molecule. The transmission results 
for [7=4 and U =8 arc depicted in the left and right panels in Fig. [3] respectively. The implications of this are two 
fold. First, both the Breit-Wigner resonance and the Fano resonance start shifting to higher energies. Meanwhile, the 
linewidth of the Breit-Wigner resonance starts decreasing whereas the Fano resonance gets broadened. We checked 
that the total linewidth stays constant. This suggests that the Breit-Wigner resonance is transfering some of its 
linewidth to the Fano resonance as U starts increasing. This is simply because the direct nonresonant electron 
transport occurring through one of the channels where the down-spin electron is located is becoming less likely as the 
Coulomb repulsion is pushing the up-spin electron away from it. Consequently, the Fano resonance, which is arising 
from the interference of the nonresonant and resonant channels, is acquiring the lost linewidth in the Breit-Wigner 
resonance. 

The transmission results in the strongly correlated limit are depicted in Fig. [4] Increasing the value of U any further 
does not alter the shape or the location of the resonances therefore we are content to show the results for only [7=128. 
It is apparent from this figure that the Fano resonance has been pushed to higher energies hence it is no longer visible. 
Moreover, the Breit-Wigner resonance in the noninteracting limit has shifted 2| (3m | to settle into the antibonding 
energy while its linewidth has been reduced to half. This behaviour is entirely different from the serially coupled case 
where the electron transport is completely inhibited due to the Coulomb blockade of the only transport channel.^ 
In parallel coupling, there are two distinct transport channels and even though the one containing the down spin 
electron is now unavailable due to Coulomb blockade, the other one still permits hopping. We believe this is precisely 
why the linewidth of the Breit-Wigner resonance has been halved compared to the noninteracting case. Moreover, the 
fact that the electron transport has shifted to the antibonding energy at this limit shouldn't come as a surprise. An 
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FIG. 4: The transmission probability T(e) for an interacting system with ?7=128 as a function of Fermi energy e. Red(solid) 
curve corresponds to Po = —10, Pm = Pcm = —2 whereas blue(dashed) curve is obtained for Pc = —10, Pm = —1 and 
Pcm = —2. 



up-spin electron hopping from left contact to the right contact uses only the atom in which the down-spin electron is 
not located without making a detour to the atom where the down-spin electron is residing. Therefore, intramolecular 
hopping is largely avoided in this regime resulting in non-overlapping wave functions for the atoms of the molecule. 
This is nothing but the antibonding configuration. 

In this paper, we investigated the electron transport through a diatomic molecular electronic device parallelly 
coupled to infinite ferromagnetic contacts at zero temperature. We modeled the system by using a model Hamiltonian 
involving Hubbard term and the contacts were described with the aid of recently developed source and sink potentials. 
Our results showed that the transmission spectrum exhibits a Fano resonance in the interacting regime as well as 
the usual Breit-Wigner resonance. We explained this as a result of the decoherence introduced by the Coulomb 
interaction with the electron residing inside the molecule. As the electron-electron interaction strength inside the 
molecule has been increased, we witnessed the gradual shift of both resonances to higher energies while the Breit- 
Wigner resonance transfered some of its linewidth to the Fano resonance. When we approached the U — > oo limit, 
only the Breit-Wigner resonance survived and got shifted by twice the coupling between the atoms in the molecule. 
We attributed this to the unavailability of the one of the transport channels due to Coulomb blockade while the other 
one remains open. This conclusion is supported by the reduction of the linewidth of the Breit-Wigner resonance to 
half of its noninteracting value. We would like to point out that the scenario that has been discussed in this paper 
can be realized experimentally. The value of the Hubbard interaction strength relative to the other couplings can be 
adjusted by stretching the junction. We expect that the U — > oo limit would be attained when the bond is about to 
dissociate. 

A.G thanks Prof. Barry Friedman for a critical reading of the manuscript and acknowledges fruitful discussions 
with Prof. Matthias Ernzcrhof. 
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